Surface-induced Magnetism Fluctuations in Single Crystal of NiBi 3 Superconductor 
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We report anistropy in superconducting and normal state of NiBig single crystals with T c = 4.06 
K. The magnetoresistance results indicate the absence of scattering usually associated with ferro- 
magnetic metals, suggesting the absence of bulk long range magnetic order below 300 K. However, 
the electron spin resonance results demonstrate that ferromagnetism fluctuations exist on the surface 
of the crystal below 150K. 
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I. INTRODUCTION 

Ferromagnetism (FM) and superconductivity (SC) are 
two fundamental condensed matter phenomena. The for- 
mer one favors all spins parallel, while in the latter one 
carriers condense in Cooper pairs where spins are an- 
tiparallcl in classical BCS scenario. The superconduct- 
ing properties of ferromagnetic superconductors are cur- 
rently under debate. The p-wave superconductivity was 
proposed to coexist with itinerant ferromagnetismii So 
far, the coexistence of FM and SC has been discovered 
in UGe2,— URhGe,— and UCoGe.— A binary intermetal- 
lic compound, ZrZn2 was reported to exhibit coexistence 
of ferromagnetism and superconductivity^ Interestingly, 
superconductivity originates from the surface alloy rich 
in Zn£ In addition, the spin-singlet pairing supercon- 
ductivity and magnetic order parameters entwine each 
other in a spatially modulated pattern, which allows for 
their mutual coexistence in for example borocarbides,— 
CeCoIn 5 in high magnetic fields,—^ HoMoeSg,— P doped 
EuFe 2 As 2r ii and ErRh 4 B 4 ^ 

Ni-Bi based compounds show complex physics phe- 
nomena including antiferromagnetism and superconduc- 
tivity. LaNiBiO is isostructural to LnFeAsO (Ln rep- 
resents the rare earth elements) iron based high tem- 
perature superconductors. CeNiBi2 is an antifcrromag- 
netic metal, whose structure is similar to that of iron 
based superconductors^ Interestingly, when Ni defi- 
cient, CcNii_ a; Bi2 shows superconductivity— Recently, 
NiBi3, an intermetallic superconductor with T c = 4.06 
K, attracted some attention, due to possible coexistence 
of ferromagnetism (Curie temperature < 750 K) and 
superconductivity in polycrystals^ Coexistence of fer- 
romagnetism and superconductivity is also observed in 
nano structured crystals below 30 KJ£ The polycrystal 
NiBi3 is reported to have large field dependent thermal 
transport properties and magneto-resistance (MR)J^ In 
addition, the magnetic ions doping of Co results in T c 
enhancement and onset of resistivity drop around 10 K 
in NiBi3.— These results suggest that magnetic order is 
related to superconductivity in NiBi3. However, other 
two groups reported that no magnetism was observed in 



NiBi3 bulk crystals . 15 ' 18 It should be noted that magnetic 
impurities cannot be avoided during synthesis process for 
nominal stoichiometric 1:3 composition according to the 
Ni-Bi phase diagram^ Thus, it is necessary to investi- 
gate the magnetism of NiBi 3 in impurity-free single crys- 
tals. Here, we report the comprehensive study of NiBi3. 
We show evidence of the ferromagnetism fluctuations in- 
duced at the surface below 160 K. No large MR or field 
dependant heat capacity is observed, which indicates that 
no bulk magnetism exists in NiBi3. 



II. EXPERIMENTAL 

Single crystals of NiBi3 were grown from self-flux 
method with Ni:Bi = 1:10 mol ratio. This method 
avoids the NiBi phase and other magnetic impurity in 
the solid state reaction with of stoichiometric Ni:Bi (1:3) 
according to the Ni-Bi phase diagram^ High purity Ni 
shot and Bi pieces were mixed and sealed in an evacu- 
ated quartz tube. The tube was heated to and soaked 
at 1150 °C for 2 hours, then cooled down to 400 °C 
with 5 °C/h. Finally, the tube was spun in a cen- 
trifuge to separate the Bi flux. As is shown in Fig. 1 
(b), needle like, silver colored single crystals with size 
of ~ 2 mm x 0.2 mm x 0.2 mm (with b axis of the 
longest dimension) were obtained. The single X-ray data 
were collected using the Bruker APEX2 software package 
Apex2 on a Bruker SMART APEX II single crystal X- 
ray diffractometer with graphite-monochromated Mo K Q 
radiation (A = 0.71073 A) at room temperature. The 
obtained lattice parameters are a — 0.8955(0.0019) nm, 
b = 0.4153(0.0012) nm and c = 1.158(0.002) nm, which 
is consistent with the previous resultsJ^ The directions 
of axis were also determined. 

Electrical resistance measurements were performed us- 
ing a four-probe configuration, with the applied current 
along the b - axis. Thin Pt wires were attached to electri- 
cal contacts made of epoxy. Electrical transport and heat 
capacity measurements were carried out in Quantum De- 
sign PPMS-9. The rotating sample holder was used to 
adjust the direction of the magnetic field. During the 
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FIG. 1. (a) Temperature dependence of pb (open circle) and 
its fitted curves (dashed line) for NiBi3. (b) The photograph 
of single crystal of NiBi3 . (c) Temperature dependence of pb 
for NiBi3 vicinity the superconducting transition. The ellipse 
marks the initial resistivity drop. 



MR measurements, H is perpendicular to the b - axis. 
Electron spin resonance (ESR) measurement was carried 
out in Bruker EMX-plus model spectrometer (v ~ 9.4 
GHz). 



III. RESULTS AND DISCUSSIONS 

Figure 1(a) shows the temperature (T) dependence of 
resistivity along the 6-axis (pb) for NiBi3 single crystal. 
Its superconducting transition can be seen in Fig. 1(c). 
The obtained T c q is about 4.06 K, which is consistent 
with former results. The superconducting transition is 
sharp with a transition width AT C ~ 0.04 K. Interest- 
ingly, the pb — T curve has a initial drop around T = 
4.3 K (Fig. 1(c)). Therefore, it rules out the possibility 
of the NiBi (S.G. PG^/mmc, NiAs type) impurity with 
T c = 4.25 K as the origin of the blunt specific heat (C) 
transition in the region from 4.1 K to 4.2 K in the poly- 
crystalline sampled The residual resistivity ratio (RRR) 
is 19.7, higher than that of the polycrystalline sampled 
The initial drop is reminiscent of filamentary supercon- 
ductivity that precedes bulk phase coherence, as observed 
in for example quasi- ID metal Nb2Se3^ The pb rises 
with positive curvature with increasing temperature be- 
low 22 K, and increases with a linear pb — T relation up 
to ~ 60 K. Then, pb increases with a negative curvature 
above 60 K, and shows a saturation tendency towards 
high temperatures. The overall Pb{T) follows the empir- 
ical model initially applied to A15 superconductors such 
as Nb 3 Sn.— The model gives 

2^ 

P = Po + Pi ■T + p 2 -cxp(— -) (1) 
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FIG. 2. (a) T 2 dependence of C/T of NiBi 3 for H = kOe 
(solid circle) and H = 20 kOe (open circle); the linear fitting 
between 3 K and 6 K for H = 20 kOe. (b)The T dependence 
of AC es = C(H = OkOe) - C(H = 20 kOe))(Left panel) and 
the T dependence of pt (right panel) for NiBi3. (c)The T 
dependence of C for MB13. The dashed line represents the 
3NR, where N and R is the number of atoms in molecule and 
universal gas constant, respectively. 



where, po is the residual resistivity; pi, P2, and To are ma- 
terial dependent parameters. The third term arises from 
the phonon-assisted scattering between the two Fermi- 
surface sheets. The fitted curve is shown in Fig. 1(a) 
as the dashed line. p = 7.453±0.016 pSl- cm, p\ = 
0.244±0.004 (til- cm K' 1 , p 2 = 96.2±1.5 /zfi- cm, and 
T = 103.8±0.7 K arc determined from the fitted results. 
These values are comparable with those of NbaSn^ 

Figure 2(a) shows the T 2 dependence of specific heat 
(C) divided by (T) of NiBi 3 single crystal for H = kOe 
and H = 20 kOe. Obviously, no magnetic specific heat 
contribution can be observed. Around 4 K, a specific heat 
jump due to superconducting transition can be observed. 
The superconducting transition can be suppressed by H 
— 20 kOe. Interestingly, if we use the Sommcrfeld - Debyc 
expression C = C e + C\ = 7T + (3T 3 (7T is the electron 
contribution, and /3T 3 is the lattice contribution) to fit 
the curve, the linear C/T — T 2 relation between 3 K and 
6 K gives unphysical negative 7. This result is consistent 
with the previous reportsJ&^i Yet, (3 = 8.255±0.02 mJ 
mol -1 K -4 is obtained from the linear fitting between 3 
K and 6 K. 9p = 98. 0±0.1K can be estimated from the 
Qd = \/12tt 4 NR/5P, where N is the number of atoms 
per molecule. In addition, the C/T — T 2 deviates from 
linearity below 3 K. Adding the anharmonic contribution 
~ T 5 docs not improve the fit. This is most likely due to 
the low Debye temperature (6d ~ 100 K)3i Most likely 
measurements of heat capacity below 2 K are necessary 
in order to reliably estimate 7. At high temperature, 
C for NiBi 3 approaches the ideal value of 3NR as the 
Dulong- Petit law predicts (Fig. 2(c)). 



3 



T = 5K 
T= 10K 
T = 20K 
T = 40K 
T = 60K 



ff- (H/p JOT)) 2 



1E-3 
10 



10° 10' 
Hp (OT) (kOe ii'cm" 1 



4.0 

3.5 

3.0 

"°2.5 
P 

S.2.0 

Q." 

a : 1.5 
< 

1.0 

0.5 

. o.o 



T = 80K 
T = 100K 
T = 150K 
T s 200K 
T = 250K 
T = 300K 



* v « • • ». . 

»v: •••• 







200 400 600 800 10 
Hip (OT) (Oe n W 1 ) 



FIG. 3. (a) MR the Kohler's log-log plot in single crystals 
of for NiBi 3 at T = 5, 10, 20, 40, 60 K. (b) MR at different 
temperatures and the Kohler's plot in single crystals of for 
NiBi 3 at T = 80, 100, 150, 200, 250, 300 K. 



In order to investigate the specific heat jump of super- 
conducting transition, C(H = kOe) - C(H = 20 kOe) is 
introduced to estimate the electronic specific heat in the 
superconducting state. Figure 2(b) shows the AC es (C(H 
= kOe) - C(H = 20 kOe)) and p b - T curves. The 
ACe S /T ~ 20 mJ moh 1 K~ 2 . From the McMillan for- 
mula 



ln( 



1.45T C 



) - 1.04 



1.04 + ln(i#2k)(l-0.62/i*)' 



(2) 



we estimate the electron-phonon coupling constant A e _ p 
~ 0.91 by assuming pf — 0.13, which is a typical value 
for the Coulomb pseudo-potential^ This indicates that 
NiBi3 is a strongly electron-phonon coupled supercon- 
ductor. 

Semiclassical transport theory predicts that Kohler's 
rule will be valid if there is a single species of charge 
carrier and the scattering time r is the same at all points 
on the Fermi surface. MR at different temperatures can 
be scaled by the expression eq.(3) with the assumption 
that scattering rate 1/r is proportional to p(T), and lo c is 
the frequency at which the H causes the charge carriers 
to sweep across the Fermi surface: 



A Pxx (H,T) 
Pxx(0,T) 



/Kr) = /( 



H 



P(0,T) 



(3) 



and the corresponding plots arc known as Kohler's plots. 

Figure 3(a) depicts the Koler's plot for NiBi3 from 5 
K to 60 K. Even at H 1 b = 9 T, the MR is only about 
4.5% at T = 5 K. This value is very typical for a nonmag- 
netic metal. All curves collapse onto a single line, sug- 
gesting that the Kohler's rule is valid (Ap xx (H) / p xx (0) 
oc (p H/p xx (0)) 2 ). Figure 3(b) depicts the MR data 
Ap xx (H )/p xx (0) as a function of poH / p xx (0) above 80 
K. Apparently, Kohler's rule is violated above 80 K. Due 
to the crystal structure, NiBi3 should be a quasi-one di- 
mensional system, whose Fermi surface topology is not 
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FIG. 4. (a)The ESR spectra of NiBi 3 measured at 1.85 K 
with sweeping H parallel to the b - axis and perpendicular to 
&-axis. Inset shows the amplified plots in the low field region. 
The ESR spectra of NiBi3 measured at different temperatures 
with sweeping H parallel to the b - axis (b) and perpendicular 
to b - axis (c). The curves are shifted for clarity. The dashed 
lines mark the maximum of the ESR spectrums. 



spherical. The possible explanation is that multiband 
conductivity in NiBi3. Below 60 K, the contribution of 
one band (with large lo c t) dominates whereas the con- 
tribution of other bands is negligible. The system can 
then be regarded as a single band conductor. Above 80 
K the contributions are comparable, leading to deviation 
of the Kohler's rule. Similar phenomenon has been ob- 
served in low dimensional system, such as LiFeP.— What 
is important, the MR is positive below 300 K, indicating 
that NiBi3 is not ferromagnetic in bulk. 

Figure 4(a) shows the ESR (dP/dH) spectra of NiBi 3 
measured at 1.85 K with sweeping H parallel to the b - 
axis (H || b) and perpendicular to 6-axis (H _L b). Ob- 
viously, anisotropy can be observed. Both spectra show 
typical ESR signals for supercondcuting state. As shown 
in the inset, a stepwise sharp signal appears at low fields 
typical for the magnetic-shielding feature, indicating the 
appearance of superconductivity. Then, ESR signal show 
a hump below H C 2 with increasing H . 

Figure 4 (b) and (c) show the ESR spectra of NiBi3 
measured at different temperatures (above T c ) with 
sweeping H \\ b and H _L b, respectively. No obvious 
resonance signal around H = 3200 Oc for paramagnetic 
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electron can be observed. Interestingly, the ESR spectra 
of NiBi 3 show abnormal diplike signal at low field below 
160 K, which can be suppressed by increasing tempera- 
ture. As marked by the dashed lines in Fig. 4 (b) and 
(c), the magnetic fields corresponding to diplike signal 
are anisotropic. In heavy fermion system CeRuPO, FM 
correlations is define by the ESR signal^ At low temper- 
ature below 10 K (Curie temperature = 15 K), ESR of 
CeRuPO shows similar diplike signals at low fields as well 
as NiBi 3 . This diplike signal should be related to ferro- 
magnetic moment fluctuations. At high temperature, the 
diplike signal of NiBi3 disappears. The microwave signal 
of ESR can only exist on the surface of a metal (in pm 
length).— This strongly suggests that ferromagnetic mo- 
ments fluctuations exist the surface on NiBi3 below 160 
K, and that NiBi3 is non-magnetic in bulk. The crystal 
structure of NiBi3 can be regarded as packing of NiBi3 
rods along the b - axis. In the bulk, the coordination of 
NiBi3 rods has perfect translational symmetry; while on 
the surface, the translational symmetry is broken. The 
FM fluctuations in NiBi 3 could be attributed to the sur- 
face effect since surface tension modifies the electronic 
band structure, favoring FM on submicron length scale 
in nanostructurcd samples 



TABLE I. Parameters for NiBi 3 . 



f(K) 




7(K) 

51.5 2.0 2.5 3.0 3.5 4.0 4.5 



E , 
O 

q 



f,J I If 



/ | f 



i i i i ' i ' i j i 

H ±c 

H=4kOe 



/ 15 ; f ' 



1 t 



(C) f 




1.5 2.0 2.5 3.0 3.5 4.0 4.51.5 2.0 2.5 3.0 3.5 4.0 4.5 

7(K) 



7(K) 



FIG. 5. Temperature dependence of the resistivity Pb{T) of 
NiBi3 for (a) H || a, (b) H \\ b and (c) H \\ c a at the various 
magnetic fields, (d) Temperature dependence of the upper 
critical field H C 2 obtained from the midpoint and 10% on the 
pt — T curves for H || a, H \\ b and H \\ c. 

Figure 5 (a), (b), and (c) show the pi, — T curves 
at the various magnetic fields for H || a, H || b 
and H || c, respectively. The upper critical field 
(H C 2) is obtained from the midpoint and 10% resistiv- 
ity on the pi, — T curves, which is depicted in Fig. 5 
(d). Obviously, H C 2 show linear temperature depen- 
dence in the experimental temperature range. The ob- 
tained [dH% id /dT]\ Tc ([dH™ % /dT}\ T J are 1.48(1.22), 
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and H \\ c, respectively. According to the conven- 
tional one-band Werthamer-Helfand-Hohenberg (WHH) 
theory^ which describes the orbital limited H C 2 of dirty 
type-II superconductors, H c2 = 0.693[dff c2 /dT] \t c T c . 
Then H^ 2 (0), H b c2 (0), ff c c 2 (0) can be obtained as 6(4.95), 
12.7(10.56), and 4.75(4.62) kOe. The data in the brack- 
ets here and the following discussion are within the 10% 
resistivity results. 

In the anisotropic superconducting system, coherence 
length (£g L , where i is the axis) can be estimated from 



H? 2 



'I'd 



• H c2 



The calcu- 



lated £§ L , £ b GL , ^ L are 18.1(18.3) nm, 38.3(39.0) nm, 
and 14.3(17.1) nm, respectively. The effective mass ten- 
sor is related to the GL coherence length as m a : : m c 

The calculated effective 



mass ratio m a : m b : m c 4.48:1:7.15(4.58:1:5.02). Ob- 
viously, NiBi3 shows a quasi one dimensional behavior 
along the b - axis, with effective mass ratio m* ~ m* and 
m*/ml ~ 5 (i represents a and c). 

The thermodynamic critical magnetic field at T 
= K (_ff c (0)) can be estimated through H c (0) — 
[47r7V(£ , F )A 2 (0)] 1/2 . According to the band calculation, 
the density of states at Fermi surface N(Ep) is 2.55 
states/eV f.u.^ By assuming the A(0) = 1.76fc B T c , H c 
is estimated to be ~ 430 Oe. This value is consistent 
with experimental rcsultsJ^ GL parameter k, the lower 
critical field at zero temperature (H c i(0)), and the pen- 
etration depth Xql at T = K can be estimated from: 



HI 



V2k*H c \H* = H c 



1iik 1 -0.18 



, and X\ 



L c2 — V ^rx, ll c ,iJ d — ±± c n GL ~ K ^GL' 

where i represents the direction of a, &, c. The estimated 



3.13(2.60), and 1.17(1.14) kOc/K for H \\ a, H \\ b, 
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parameters in the superconducting state and the param- 
eters in the normal state mentioned above are listed in 
Tab. I. 



IV. SUMMARY 

We investigated anisotropy in superconducting and 
normal state properties of N1EU3. The carrier scattering 
mechanism is dominated by electron-phonon scattering 
whereas MR suggests the absence of magnetic scatter- 
ing commonly observed in bulk FM materials. The heat 
capacity results indicates it is a strongly e — p coupling 
superconductor. However, FM fluctuations are detected 
at the surface of the crystal. Finally, we give anisotropic 



superconducting parameters. 
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